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Preface 


It  is  suggested  that  debonding  and  fiber  fracture  in  composite 
materials  might  lead  to  significant  internal  energy  losses.  In  the 
past,  this  damage  has  been  assessed  primarily  by  nondestructive 
techniques  such  as  ultrasonics  and  radiography.  This  was  an 
experimental  study  to  test  the  feasibility  of  using  the  measurement  of 
internal  energy  loss  as  a  nondestructive  technique  to  determine  the 
amount  of  damage  that  has  occurred  in  a  composite  during  a  fatigue 
test.  I  performed  fatigue  tests  on  graphite-epoxy  composite  specimens 
in  order  to  quantify  the  damage  in  terms  of  the  damping  characteristics 
of  the  composite. 
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through  radiography  and  ultrasonics. 
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support.  I  am  deeply  grateful  to  Dr.  Peter  Torvik,  my  thesis  advisor, 
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understanding  throughout  this  study. 
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Abstract 

The  purpose  of  this  experimental  study  was  to  investigate  the 
possibility  of  quantifying  damage  in  composites  in  terms  of  damping 
characteristics.  Bending  fatigue  tests  were  performed  on 

((0/-45/90)s)2  T300/5208  graphite-epoxy  composite  specimens.  Testing 
was  performed  by  exciting  a  cantilever  beam  so  as  to  maintain  a  con¬ 
stant  amplitude  at  resonant  frequency.  Results  indicate  that  the 
amount  of  damage  occurring  in  this  composite  can  be  explained  in 
terms  of  an  increase  in  its  material  damping.  The  resonant  fre¬ 
quency  was  found  to  decrease  as  the  material  damage  increased. 
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MATERIAL  DAMPING  AS  A  MEANS  OF 
QUANTIFYING  FATIGUE  DAMAGE  IN  COMPOSITES 

I.  Introduction 


Background 

Composite  materials  are  ideal  for  structural  applications  where  high 
strength-to-weight  and  stiffness-to-weight  ratios  are  required.  Aircraft 
and  spacecraft  are  typical  weight-sensitive  structures  in  which  composite 
materials  are  cost-effective. 

Composites  designated  as  advanced  fiber-reinforced  composite 
materials  are  those  with  new  ultrahigh  strength  and  stiffness  fibers 
such  as  boron  and  graphite.  The  matrix  material  can  be  a  plastic 
such  as  epoxy  or  a  metal  such  as  aluminum.  Advanced  composites  have 
two  major  advantages:  improved  strength  and  stiffness. 

Composites  can  be  made  three  times  as  strong  as  aluminum  yet  weigh 
only  60  percent  as  much.  Other  composites  can  be  made  that  have  the 
same  strength  and  stiffness  as  high-strength  steel,  but  are  70  percent 
lighter.  Moreover,  composite  materials  can  be  tailored  to  efficiently 
meet  design  requirements  of  strength,  stiffness,  and  other  parameters. 
These  advantages  should  lead  to  new  and  more  efficient  aircraft  and 
spacecraft  designs  in  contrast  to  those  based  on  conventional  materials. 

The  strength  of  composites  is  governed  by  their  flaw-initiated 
characteristics.  Thus,  the  mechanics  of  fracture  including  crack 
propagation  or  extension  are  of  extreme  importance  in  the  design 
analysis  of  composite  structures.  The  fracture  process  takes  place 
in  three  stages.  First,  an  initial  flaw  may  be  present  or  a  micro¬ 
crack  initiated.  Second,  the  microcrack  grows  and  may  link  with 
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others  to  attain  microcrack  size.  Third,  at  a  critical  stress  level,  the 
crack  propagates  through  the  material.  Fracture  occurs  when  the  crack 
reaches  a  critical  length. 

Salkind  (Ref  13)  has  proposed  that  composites  may  fail  by  matrix 
cracking,  fiber  failure,  and  delamination,  either  separately  or  in 
combination.  He  also  suggests  that  a  fatigue  failure  criterion  could 
be  the  number  of  cycles  to  a  given  change  in  stiffness  rather  than 
the  number  of  cycles  to  fracture  since  the  loss  of  stiffness  can  result 
in  structural  failure  long  before  complete  fracture  occurs. 

The  internal  damping  of  fiber-reinforced  composites  is  more 
sensitive  than  stiffness  to  the  failure  mechanisms  previously  mentioned. 
Accordingly,  damping  measurements  can  be  valuable  in  the  early  detection 
of  structural  damage. 

Much  of  the  experimental  work  on  the  damping  and  stiffness 
properties  of  composite  materials  has  been  concerned  with  measuring 
frequency  and  temperature  dependence  (Ref  4) .  Both  aluminum  and  fiber- 
reinforced  plastic  beams  have  been  tested  by  the  methods  of  logarithmic 
decrement  and  resonant  dwell  (Ref  8) .  Air  damping  was  shown  to  be 
significant  at  low  frequencies  and  large  amplitudes.  Schultz  and 
Tsai  (Ref  14)  used  logarithmic  decrement  and  bandwidth  techniques  to 
test  aluminum  and  glass  fiber-reinforced  plastics  in  flexural 
vibration.  They  found  that  the  amount  of  fatigue  cracking  correlated 
reasonably  well  with  the  amount  of  increase  in  damping.  Ulitchny  and 
Lucas  (Ref  10)  studied  the  change  in  the  damping  of  free  vibrations 
and  found  crack  damage  present  as  changes  in  damping  occurred. 

For  the  commercial  utilization  of  composites,  it  is  necessary 
to  ensure  material  integrity  through  reliable  methods.  Damage 
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evaluation  by  nondestructive  techniques  (NDT)  is  currently  being 
investigated.  Of  the  NDT  methods  available,  holography  and 
ultrasonics  show  the  most  promise.  Because  damping  can  be  expected  to 
be  sensitive  to  microstructural  damage,  damping  measurements  may  provide 
a  means  of  achieving  the  necessary  quality  control. 

This  thesis  studied  the  feasibility  of  using  the  change  in 
material  damping  to  quantify  structural  damage  in  the  course  of  a 
destructive  bending  fatigue  test  on  composite  beams.  Observed 
changes  in  damping  were  compared  with  NDT  results  in  order  to  provide 
a  comparison  between  their  two  means  of  evaluating  damage. 


II.  Analysis  of  Beam  Stress  and  Displacement 


Beam  Theory 

The  composite  specimens  were  modeled  as  beams  and,  as  such,  a 
few  simplifying  assumptions  have  been  made: 

1)  Beam  deformations  are  described  by  the  elementary  theory 

of  flexure,  so  that  plane  sections  normal  to  beam  axis  remain 
plane  and  normal  to  the  deformed  axis 

2)  Each  ply  is  assumed  to  satisfy  the  linear,  elastic  constitutive 
relationship 

3)  Flies  are  assumed  to  have  uniform  thickness 

The  well-known  differential  equation  of  t hi  !  r  <  >  ->  urve  of  a 
beam  whose  vibrations  occur  in  its  plane  of  svwn-t  •  .  f  •> )  is 

(1) 

ax 

in  which 

y (x)  =  deflection  of  the  neutral  axis, 

El  =  stiffness,  and 

M  =  bending  moment  at  any  cross-section.  The  direction 
of  the  axis  and  positive  directions  of  bending  moments  and  shear 
forces  are  shown  in  Figure  1. 


Differentiating  Eq  (1)  twice  gives 


(4) 


(4a) 
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with  time  and  can  be  expressed  as 

y  »  X(A  cospt  +  B  sinpt)  (5) 

Where  X  is  a  function  of  the  coordinate  x.  Substituting  Eq  (5)  into 
Eq  (4)  leads  to 


C  C  C  C 

where  1,  2,  3,  and  4  are  constants  to  be  determined  from  the  end 

conditions  of  the  beam. 

Assuming  that  one  end  (x  =  0)  is  built  in  and  the  other  end  is 
free,  the  following  end  conditions  are  obtained: 
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The  other  two  give 


The  first  two  conditions  result  in  Cj*  0  . 

the  following  equation  from  which  the  frequency  can  be  obtained: 

(coskl) (coshkl)  *  -1 

When  a  tip  mass  is  added  to  the  previously  free  end  of  the 
beam,  new  end  conditions  are  obtained: 


Where  m  =  tip  mass.  The  rotatory  inertia  of  the  tip  mass  has  been 
neglected.  Again,  the  first  two  conditions  give  cxsr  C3*  0  in 
Eq  (8).  The  last  two  give  the  following  equation: 

(coskl) (coshkl)  +  1  p2m  ^ 

(sinkl) (coshkl)  -  (coskl) (sinhkl)  k^EI 

From  Eqs  (7)  and  (4a) 

2  EIkA  (9a) 

p  “  P  A 

and  Eq  (9)  can  be  expressed  as 
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(coskl) (coshkl)  +  1 


(sinkl) (coshkl)  -  (coskl) (sinhkl) 


(9b) 


This  transcendental  equation,  together  with  Eq  (9a),  must  be  solved 
In  order  to  find  the  natural  frequencies  in  terms  of  the  mass  ratio 
of  tip  mass  to  beam. 

An  approximate  method  known  as  Rayleigh's  Quotient  can  also  be 
used  to  estimate  the  frequency  of  vibration.  Its  application  comes 
from  consideration  of  the  total  energy  of  the  system.  In  applying 
this  method,  an  initial  assumption  has  to  be  made  regarding  the 
beam  deflection  shape  during  vibration.  In  this  case,  the  deflection 
was  taken  to  be 

X  -  f(x)  -  x2(31-x)  (10) 

The  potential  energy  of  the  beam  is  simply  its  strain  energy: 

/  l  ^(S) 2  (sinpt)2 

As  the  tip  mass  is  assumed  to  contribute  no  potential  energy,  the 
maximum  potential  energy  is  then 

The  kinetic  energy  in  the  beam  is 
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('5t)  dx  "  J />Ap2X2cos2pt  dx 
This  gives  a  maximum  kinetic  energy  of 


for  the  beam.  The  tip  mass  contributes  to  the  kinetic  energy  of 
the  system: 


Then  the  maximum  kinetic  energy  of  the  tip  mass  is 


m  2  2 
2  P  * 


X“1 


(13) 


The  total  energy  of  the  system  is  a  constant  and  equal  to  the 
sums  of  the  kinetic  and  potential  energies.  It  follows  that  the 
maximum  kinetic  and  potential  energies  of  the  system  can  be 
equated.  This  yields: 


f  ei  / dzx\2d  r 

/  "~2  IT?/  ^  / 

J%  \dx  /  Ja 


PA  2  2, 
-‘j-  p  X  dx 


4.  m  n2v2 

+  2  P  X 


x=l 


Rearranging,  the  frequency  can  be  found  to  be: 


(14) 


Substituting  Eq  (10)  into  Eq  (14)  gives  the  frequency  in  terms  of 
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the  assumed  beam  deflection: 


(61- 6x)  dx 


x4(31-x)2dx  +  [f  x2(31-x)2Jx-1 


carrying  out  the  integration  yields 


2 

P 


33  pAl 


210  El 

. 


In  both  Rayleigh's  Quotient  and  the  exact  solution,  p  can  be  solved 

for  in  terms  of  the  tip  and  beam  masses.  The  necessary  material  constants 

for  the  material  used  in  this  investigation  were  determined  and  are  listed 

in  Table  I.  Eq  (9)  can  be  solved  for  kl  by  an  iterative  process  and  the 

resonant  frequency,  V L,  can  be  determined  from  Eq  (9a).  Thi^  yields 

K 

p  =  174.96  rad/sec  or  V  **  27.85  cps  for  exact  solution.  Using  Rayleigh's 

K. 

method,  p=171.62  rad/sec  or  *  27.31  cps.  There  is  less  than  two  per¬ 
cent  error  in  the  frequency  determined  by  both  methods. 

Root  Stress  and  Tip  Amplitude 

Since  the  experimental  procedure  employed  in  this  investigation  used 
the  measurement  of  tip  displacement  as  a  means  of  controlling  beam  stress, 
it  was  necessary  to  develop  a  relationship  between  stresses  and  tip 
displacements.  Since  the  beam  was  assumed  to  be  in  simple  bending,  nominal 
beam  stresses  were  taken  to  be 


<r.«£ 

i 


Incorporating  Eqs  (1)  and  (5),  the  maximum  outer  fiber  stress  at 
any  point  along  the  beam  is 
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Table  I 


i 


Specimen  Material  Constants 


Young's  Modulus,  E^ 

Area  Moment  of  Inertia,  I 
Stiffness,  El 


1.138  x  107  psi 
5.1177  x  10-5  In4 
5.8233  x  102  lb-in 


Beam  thickness,  2c 
Beam  width,  D 

Effective  beam  length  from 
root  to  tip  mass,  1 

Beam  Mass,  pAl 

Tip  Mass,  m^,  includes  accelero¬ 
meter,  fiberglass  doublers, 
beam  tip  beyond  center  of 
accelerometer 

Mass  Ratio,  “T 
pAl 


0.085  in 
1.00  in 
5.12  in 

7.734  x  10  4  slugs 
4.446  x  10  3  slugs 


5.749 


1.  Young's  Modulus  was  computed  by  the  method  described 
in  Appendix  A 
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2 

O  *  cE  d  X  (Acospt  +  Bsinpt) 
dx2 

The  stress  at  the  root  (x=*0)  is  then 
2 

=  cE  d  X(0)  (Acospt  +  Bsinpt) 

j  2 

dx 

Acceleration  of  the  tip  can  be  expressed  as 
d2  2 

a  =  — 2.  =  p  x(Acospt  +  Bsinpt) 

T  dt2  x=l 


(16) 


(17) 


Since  the  ratio  of  root  curvature  to  tip  displacement  is 
determined  solely  by  the  mode  shape,  it  is  convenient  to  combine  Eqs  (16) 
and  (17)  into 


m  £l  *r 

p2  X(l) 


(18) 


Using  the  approximate  mode  shape,  Eq  (10),  we  find 


0„  m  cE  f"(0) 

p2  f(» 


(19) 


Since  the  tip  amplitude  is  related  to  the  acceleration  through 


a 


T 


(20) 


the  tip  amplitude  required  to  produce  a  desired  stress  at  the  root 
i8 


aT  _  g(D 

cEf " (0) 


cE  3 


(21) 
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Using  the  exact  mode  shape,  Eq  (8),  with  =  0  and 


C2  "  ~C4(sinkl  +  sinhkl) 
coskl  +  coshkl 


(22) 


results  in 


X(x)  =  C4 


(sinkx-sinhkl) 


(sinkl  4-  sinhkl) (coskx-coshkx) 
coskl  +  coshkl 


(23) 


Thus 


X(l)  =  2C4 


(sinkl) (coshkl) -(coskl) (sinhkl) 
coskl  +  coshkl 


(24) 


and 


X"(0) 


2C4k2 


sinkl  +  sinhkl 
coskl  +  coshkl 


Thus 


a 

R 


ceM  k2  J"  sinkl  +  sinhkl _ 

p2  L  (sinkl) (coshkl)- (coskl) (sinhkl) 


For  the  beams  used  in  this  investigation,  Eq  (9b)  was  found  to 
yield  a  value  of  kl  =  .84142  or 


-  F1 

cE  L 


sinkl) (eoshkl)-(coskl) (sinhkl) 
2  (sinkl  +  sinhkl) 


)" 


or 

aT  -  .3312  ffRl2 
cE 
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The  formula  based  on  the  approximate  mode  shape  differs  from  this 
by  less  than  one  percent. 

For  a  beam  having  the  properties  given  in  Table  I, 


<fR  =  55,713  6t  psi/in 


Inclusion  of  a  hole  centered  a  distance  b  from  the  root  leads 
to  stresses  at  the  hole  edge  (Ref  12)  given  by 


®max  =  f  6M(b) 

c2(D-2r) 


Where  f  =  1.79  + 


1  +  — 
c 


1  +(4  1  +(rf  J 

1  -  1.04  +  1.22  (~j  J 


r  =  hole  radius 


D  =  width  of  beam 


2c  =  thickness  of  beam 


M  *  moment  at  hole  =  Ely" 


The  stress  at  the  hole  can  be  expressed  in  terms  of  the 
nominal  stress  where  the  nominal  stress  is 


ff  =  M(b)c 
n 


6M(b) 
„  2 


The  stress  concentration  factor  is  found  to  be 
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and  has  a  value  of  2.0701  for  r=0.098  in,  D=1.00  in,  and  c=0.085  in. 
The  nominal  stress  can  be  written  in  terms  of  the  tip  amplitude 

since 


M  =  EIX*'  (hole) 


so  that 


M(b)  =  El 


d2X(b) 

dx2 


Substituting  the  exact  mode  shape  given  by  Eq  (23)  and  eliminating 
through  the  use  of  Eq  (24),  the  nominal  stress  is  found  to  be 


Ec  dT 


(sinkb-sinhkb) (coskl+coshkl)-  (sinkH-sinhkl) 
2  [(sinkl) (coshkl)- (coski) (sinhkl)] 


(sinkl+sinhkl) (coskb-coshkb) _ 

2  [(sinkl) (coshkl)- (coskl) (sinhkl)] 


For  a  beam  having  the  properties  in  Table  I, 


a  *  25,089  6~  psi/in 
n  T 


Combining  this  with  the  stress  concentration  factor  given  by 
Eq  (27),  the  stress  at  the  edge  of  the  hole  is  found  to  be 

amax  =  52,384  psi/in 

where  <5^  is  the  tip  displacement.  This  value  is  seen  to  be  slightly 
less  than  the  stress  at  the  root  of  the  beam;  however,  failure  may  be 
expected  to  initiate  at  the  hole  edge  because  of  inter-laminar  edge 
stresses  (Ref  14)  and  imperfections  generated  by  the  process  of 
creating  the  hole. 
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III.  Experimental  Procedures 

Specimen  Description 

Twelve  specimens  with  the  geometry  shown  in  Figure  2  were 
supplied  by  the  Air  Force  Flight  Dynamics  Laboratory.  They  were 
machined  from  a  single  16-ply,  Thornel  300/5208,  graphite-epoxy  panel 
fabricated  by  Monsanto  Research  Corporation.  The  multiple  layers  are 
arranged  symmetrically  about  the  middle  surface  at  orientations  of 
0,  45,  -45,  90,  90,  -45,  45,  0  degrees  to  the  x-axis  starting  with 
the  top  ply  and  progressing  to  the  middle  surface.  These  orientations 
are  repeated  for  the  remaining  eight  plies  below  the  middle  surface. 

This  layup  can  be  expressed  as  ((0/-^45/90  )  )2  .  Rectangular 
fiberglass  doublers,  .118  inch  thick,  with  one  edge  beveled  at  30 
degrees  were  adhesively  bonded  to  the  specimens.  A  circular  hole  was 
drilled  in  the  center  of  each  specimen.  Typical  dimensions  of  the 
specimens  were  one  inch  wide,  0.085  inch  thick,  eight  inches  long, 
with  a  center  hole  diameter  of  0.20  inch.  An  aluminum  beam  of  the 
same  geometry  as  the  specimens  was  used  to  develop  the  testing  pro¬ 
cedures  as  was  an  unusable  specimen,  AB  70.  After  mounting,  the 
effective  length  of  each  specimen  from  root  to  tip  was  5.12  inches  as 
defined  in  Figure  2. 

Constant  Amplitude  Cycles 

Each  specimen  was  taken  to  failure  where  failure  for  this 

experiment  was  defined  as  extensive  visual  fiber  delamination  accompanied 

by  a  significant  change  in  resonant  frequency.  Stress  levels  were 

4  6 

chosen  so  that  failure  would  occur  between  10  and  10  cycles.  From 
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these  stresses,  corresponding  tip  amplitudes  were  calculated  from 
Eq  (21).  Constant  tip  amplitude  was  maintained  in  order  to  keep  a 
constant  maximum  stress  at  the  root  of  the  beam  while  at  resonance. 
Resonance  was  found  by  observing  the  peaking  of  beam  response  with 
varying  excitation  frequency  on  an  oscilloscope. 

An  M-B  electrodynamic  vibration  exciter,  model  CIOE,  was  used 
for  all  tests.  As  the  schematic  in  Figure  3  shows,  it  was  driven 
by  an  M-B  automatic  vibration  exciter  control,  model  N695/696,  via 
an  M-B  power  amplifier.  The  moving  element  of  the  exciter  was 
fitted  with  an  aluminum  fixture  that  served  as  the  mount  for  each 
specimen.  Figure  4  illustrates  this  arrangement.  An  Endevco  piezo¬ 
electric  accelerometer,  model  2235C,  was  glued  onto  the  specimen 
free  end  and  provided  feedback  to  the  control  unit  through  two 
Endevco  preamplifiers,  models  2614B  and  2621.  The  control  unit 
integrated  the  acceleration  twice  and  drove  the  shaker  with  a  sinu¬ 
soidal  signal  so  as  to  shake  the  specimens  vertically  and  maintain 
the  specimens  at  constant  tip  amplitude  at  resonance.  This  ampli¬ 
tude  was  read  and  monitored  directly  from  the  control  unit.  The 
same  accelerometer  and  preamplifiers  sent  a  signal  to  a  Tektronix 
oscilloscope,  model  465M,  which  monitored  the  free  end  acceleration. 
Another  accelerometer  of  the  same  type  was  mounted  on  the  aluminum 
block  at  the  clamped  end  of  the  specimens.  Its  signal  was  fed  into 
an  M-B  vibration  meter,  model  N524,  from  which  base  displacement 
was  monitored. 

A  CMC  universal  counter-timer,  model  726C,  monitored  the 
driving  frequency  continuously.  Cycles  were  counted  using  this 
resonant  frequency  for  each  beam.  This  frequency  and  the  time  of  each 


Figure  3*  Sohematio  of  Vibration  Test  System 


test  run  were  used  to  determine  the  number  of  cycles  applied  in  each 
specimen. 

Twisting  was  not  pronounced  enough  to  be  noticed  in  the  system. 
Logarithmic  Decrement 

In  a  system  with  ^<1,  the  degree  of  damping  can  be  defined  in 
terms  of  successive  peak  values  of  a  tip  amplitude  as  a  function  of 
time  obtained  in  free  vibration.  For  any  two  maxima  separated  by 
n  cycles  of  oscillation  as  shown  in  Figure  5,  the  ratio  of  these 
two  maxima  is 


=  e 


-2nT]n/(l-  T)2)** 


from  which  rj  can  be  expressed  as 


x 

r 
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-2  7Tn 

A  Honeywell  visicorder,  model  1508,  using  a  Honeywell  galvonometer, 
M-40-350A,  recorded  the  free  vibration  response  of  the  specimen  tip 
each  time  a  damping  measurement  was  made  by  logarithmic  decrement.  This 
response  signal  was  provided  to  the  visicorder  by  the  tip  accelerometer 
through  the  preamplifiers.  The  damping  measurements  using  logarithmic 
decrement  were  all  taken  at  a  root  stress  level  of  <7^  *  16450  psi. 
Bandwidth 


In  the  bandwidth  method,  damping  is  indicated  by  the  width  of  the 
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response  curve  near  the  resonant  frequency,  v  ,  as  shown  in  Figure  6. 

K 

Designating  the  width  as  a  frequency  increment,  measured  at  the 

half-power  points  (at  a  value  of  R=^max/ >/2) >  the  damping  of  the 

system  is  defined  as  =  for  values  of  ^<0.1, 

1  R 

The  tip  accelerometer  through  the  preamplifiers  provided  a 
signal  to  the  oscilloscope  to  monitor  the  free  end  response  of 
each  specimen.  This  response  was  recorded  in  order  to  determine 
damping  by  the  bandwidth  method  as  frequency  varied  with  constant 
base  displacement. 

Tables  III  through  XIII  illustrate  the  intervals  at  which  damping 
measurements  were  made  for  each  specimen  using  both  bandwidth  and 
logarithmic  decrement.  Initially,  logarithmic  decrement  only  was 
used  to  calculate  damping.  The  bandwidth  technique  was  not  incorpor¬ 
ated  until  the  reading  at  10^  cycles  since  the  specimens  had  to  be 
vibrated  during  the  bandwidth  measurements.  At  10^  cycles  and  subse¬ 
quent  intervals  of  measurement,  the  number  of  cycles  put  on  the  beam 
during  the  measurement  were  negligible  compared  to  the  total  number 
of  cycles  on  the  beam.  Bandwidth  damping  measurements  were  taken 
at  a  representative  root  stress  level  which  occurred  at  the  stress 
level  midway  between  the  root  stress  at  the  peak  response  and  the 
root  stress  at  the  half-power  points  of  the  response  curve. 

NDT 

The  Air  Force  Materials  Laboratory  performed  nondestructive 
tests  on  the  specimens  after  the  fatigue  tests  to  determine  the 
extent  of  damage.  Radiographic  inspection  was  performed  using  a 
Norelco  150  KVP  unit.  The  specimens  were  treated  with  tetrabrome- 
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thane  on  both  front  and  back  sides  for  thirty  minutes  prior  to 
conducting  the  x-ray  inspection.  X-rays  were  taken  at  0,  —  15 
degrees  from  the  perpendicular. 

Ultrasonic  C-scans  were  taken  by  a  Sperry  UM  721  Ref lectoscope. 
The  ultrasonic  frequency  range  of  the  instrument  amplifier  was 
ten  to  twenty  megahertz.  The  ultrasonic  attenuation  of  the  specimens 
was  four  decibels  with  respect  to  water. 
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IV.  Results  and  Discussion 


Bending  Fatigue 

Results  of  the  bending  fatigue  tests  are  presented  in  Table  II. 
The  stress,  <*max,  is  the  stress  computed  to  exist  at  the  edge  of  the 
hole  from  the  observed  tip  displacement.  The  fatigue  data  in  terms 

of  maximum  root  stress  have  been  plotted  in  the  form  of  an  S-N  curve 

/ 

in  Figure  7.  The  line  drawn  through  the  experimental  points 
represents  the  best  fit  by  visual  means.  Tests  on  specimens  AB  96 
and  AB  98  were  terminated  due  to  mechanical  failure.  AB  70  was  not 
tested  due  to  the  extensive  damage  present  prior  to  the  experiment. 

Single,  continuous  logitudinal  cracks  were  observed  to  initiate 
at  the  hole  and  progress  to  the  root  of  each  specimen.  The  width  of 
the  damage  zone  was  equal  to  the  diameter  of  the  hole.  Figures  8  and 
9  illustrate  damage  in  three  of  the  beams  which  was  typical  of  the 
damage  in  each  of  the  specimens.  Although  not  presented  here,  the 
x-rays  taken  at  *15  degrees  from  the  perpendicular  indicated  that 
the  damage  was  not  limited  to  the  outer  layer. 

Damping 

Material  damping  data  were  taken  at  selected  cycle  increments 
and  calculated  from  logarithmic  decrement  and  bandwidth.  As  the 
number  of  cycles  on  each  beam  increased,  the  damping  increased  as 
expected  due  to  the  permanent  microstructural  damage  occurring  in  the 
specimens.  The  damping  results  show  the  same  general  pattern  of 
Increase  for  each  specimen  regardless  of  root  stress  level  and  total 
number  of  cycles  applied.  Figure  10  typically  illustrates  how  the 
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Table  II 


Fatigue  Test  Results 

Specimen 

Maximum  Stress  at  Hole, 
‘'max’  ?sl 

Root  Stress, 
aR,  psi 

Cycles  to 
Failure 

AB  74 

47850 

50700 

1.07  x  106* 

AB  77 

51530 

54600 

6.83  x  105 

AB  81 

63500 

67410 

9.63  x  104 

AB  82 

70120 

74660 

3.96  x  104 

AB  90 

74110 

78560 

3.01  x  105 

AB  93 

64430 

68530 

2.39  x  105 

AB  94 

58150 

61840 

6.93  x  105 

AB  95 

75960 

80780 

1.42  x  105 

AB  96** 

53430 

56830 

AB  98** 

53960 

57380 

AB  99 

51860 

55160 

2.80  x  106 

*  Failure  not  observed;  test  terminated 

**  Test  terminated  due  to  mechanical  failure 
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Log  Base  10  of  Cycles  to  Failure 


Figure  9 


Ultrasonic  C-scans 
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Figure  10 


damping  changed  as  the  number  of  cycles  increased.  Each  damping  value 
plotted  was  calculated  by  logarithmic  decrement  at  the  same  stress  level, 
<^=16450  psi.  Logarithmic  decrement  and  bandwidth  damping  values  for 
all  the  beams  are  tabulated  in  Tables  III  through  XIII.  The  values  in 
parentheses  next  to  bandwidth  in  these  tables  are  root  stresses  in  psi 
at  which  the  bandwidth  measurements  were  made. 

An  experiment  was  performed  in  order  to  determine  how  much  of  the 
total  damping  was  material  damping  and  how  much  was  air  damping.  As 
shown  in  Figure  11,  an  aluminum  paddle  was  attached  to  the  free  end 
of  one  of  the  specimens,  centered  on  the  beam  and  flush  with  the  tip; 
the  total  surface  area  was  increased  by  a  factor  of  two.  Free 
vibration  measurements  were  taken  of  the  beam  by  itself  and  then  of 
the  beam  with  the  paddle.  Air  damping  was  determined  as  explained  in 
Appendix  C.  Results  show  that  air  damping  was  present  but  was  only  on  the 
order  of  five  percent  of  the  total  damping.  Therefore,  the  damping 
results  are  primarily  indicative  of  material  damping  behavior. 

In  general,  the  bandwidth  and  logarithmic  decrement  damping 
values  were  in  fair  agreement.  This  may  be  due  to  taking  damping 
measurements  at  different  stress  levels;  bandwidth  measurements  were 
higher  and  taken  at  higher  stress  levels  than  those  measured  by 
logarithmic  decrement. 

For  each  logarithmic  decrement  curve  obtained,  damping  values 
varied  as  the  stress  decreased  while  the  beam  damped  out  in  free 
vibration. 

Damping  calculated  by  logarithmic  decrement  increased  as  a 
greater  root  stress  was  applied  at  the  same  number  of  cycles. 
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Table 

III 

Specimen  AB74 

0^  =  50700  psi 

Rel  Hum  =55%  T  -  74°F 

N 

Log  Dec 

Bandwidth (367  30) 

0 

.0035 

105 

.00547 

106 

.00493 

.00996 

Table 

IV 

Specimen  AB77 

o  =  54600  psi 

K. 

Rel  Hum  =52%  T  =  74°F 

N 

Log  Dec 

Bandwidth (50170) 

0 

.00424 

10* 

.00494 

4 

3  x  10 

.00494 

io5 

.00563 

.00677 

3  x  105 

.00602 

.00714 

6.83  x  105 

.00593 

.00604 

8.73  x  105 

.00627 
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Table  V 


) 


1 

1 

I 


f  i 

i  ! 


' 

Specimen  AB81 

=»  67410  psi 

Rel  Hum  =52%  T  ■  74°F 

► 

.  : 

N 

Log  Dec 

Bandwidth (56330) 

t' 

0 

.00420 

%-  ' 

»r 

r 

4 

1° 

.00426 

m 

4 

3  x  10 

.00565 

r 

Jt 

9.63  x  104 

.00541 

.00660 

t 

Table  VI 

Specimen  AB82 

O  =  74660  psi 

K 

Rel  Hum  =54%  T  =  67°F 

■>- 

;  n 

Log  Dec 

Bandwidth (64250) 

,  o 

.00270 

104 

.00561 

3.96  x  104 

- 

.00605 

.00471 

1 

i 
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Table  VII 


Specimen  AB90 

<7r  =  78560  psl 

Rel  Hum  -  50%  T  *  63°F 

N 

Log  Dec 

Bandwid  th(63370) 

0 

.00379 

10A 

.00542 

3  x  104 

.00553 

2.52  x  105 

.00735 

3.01  x  105 

.00710 

.00823 

Table  VIII 

Specimen  AB93 

a  =  68530  psl 

K 

Rel  Hum  =45%  T  =  82°F 

N 

Log  Dec 

Bandwidth (57210) 

0 

.00350 

io4 

.00509 

3  x  104 

.00515 

105 

.00511 

2.39  x  105 

.00602 

.00567 

36 
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Table  XI 


Specimen  AB96 

a  =  6830  psi 

R 

Rel  Hum  *  44% 

T  =  75°F 

N 

Log  Dec 

0 

.00441 

io4 

.00591 

3  x  104 

.00587 

. 

Table  XII 

Specimen  AB98 

57380  psi 

R el  Hum  =  41% 

T  =  74°F 

N 

Log  Dec 

0 

.00443 

104 

.00516 

3  x  104 

.00523 

105 

.00526 

00526 


Table  XIII 


Specimen  AB99 

o  =  51860  psi 

K 

Rel  Hum  =40 Z  T  *  81°F 

N 

Log  Dec 

Bandwidth (42250) 

0 

.00468 

V 

.00510 

3  x  104 

.00559 

io5 

.00590 

.00614 

3  x  105 

.00580 

.00688 

io6 

.00596 

.00522 

2  x  106 

.00572 

.00522 

2.80  x  106 

.00580 

.00522 
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V 


Figure  12  exhibits  typical  results  of  the  damping  increase  with  increase 
in  root  stress.  The  same  results  occurred  for  each  cycle  increment  at 
which  damping  was  measured. 

Frequency 

For  each  specimen,  the  initial  resonant  frequency  dropped  to  a 

3 

lower  resonant  frequency  within  3  x  10  cycles.  Resonance  occurred 
at  succeedingly  lower  frequencies  throughout  the  test  runs.  Typically, 
the  final  resonant  frequency  was  no  more  than  three  cycles  less  than 
the  initial  value  as  shown  in  Figure  13  and  Table  XIV.  This  decrease 
is  due,  as  is  the  damping  increase,  to  the  permanent  damage  experienced 
by  the  specimens.  Since  frequency  is  a  measure  of  stiffness,  this 
reduction  indicates  a  permanent  change  in  the  stiffness  of  the  material. 
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Log  Base  10  of  Cycles  to  Failure 
Figure  13.  i^-N  Plot,  Specimen  AB  94 


Table  XIV 


Natural  Frequency,  CPS 


Specimen 

0  Cycles 

4 

10  Cycles 

10^  Cycles 

10^  Cycles 

At  Failure 

AB  74 

26.1 

25.6 

25.6 

* 

AB  77 

27.5 

26.8 

26.6 

25.9 

AB  81 

28.2 

27.3 

26.6 

AB  82 

27.9 

26.7 

25.6 

AB  90 

26.6 

25.9 

25.6 

24.9 

AB  93 

27.0 

26.2 

25.8 

24.7 

AB  94 

28.6 

27.7 

27.6 

26.6 

AB  95 

27.8 

26.8 

25.6 

25.0 

AB  96 

28.3 

27.4 

** 

AB  98 

26.9 

26.1 

26.0 

** 

AB  99 

28.5 

27.9 

27.7 

26.8 

26.8 

* 

/ 


*  Test  terminated;  no  failure  observed 
**  Test  terminated  due  to  mechanical  failure 
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V.  Conclusions 


It  is  concluded  that: 

1)  Test  results  indicate  that  it  is  feasible  to  quantify 
damage  in  composites  in  terms  of  changes  in  material 
damping;  material  damping  increased  as  damage  occurred 
in  the  graphite-epoxy  specimens. 

2)  Frequency  decreases  as  damage  increases.  This  a  measure 
of  the  material  stiffness  which  may  also  prove  to  be  a 
valid  measure  of  damage. 

3)  Measured  material  damping  is  dependent  on  stress  level; 
nigher  stresses  yield  large  damping  values. 

4)  Air  damping,  though  present,  was  not  a  significant  por¬ 
tion  of  the  measured  damping.  Therefore,  the  damping 
values  reflect  a  true  change  in  the  material  damping. 
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Appendix  A 

Development  of  the  Composite  Stiffness 

In  order  to  determine  the  stiffness  of  the  composite,  it  is 
necessary  to  use  the  moment  equation: 

M  m  f  O  z  dx 
x  J_c  x 

where  is  the  moment  per  unit  length  as  shown  in  Figure  A-l  and 

O ^  can  be  expressed  in  terms  of  the  laminate  middle  surface  strains, 

€°,  and  curvatures,  K 
x’  x 

[«„  ]„{«;}  +  *K) 

Then  the  moment  can  be  written  in  terms  of  this  stress-strain 
relationship: 


where  Q  are  transformed  reduced  stiffnesses  and  z^>  z^  ^  are  defined 

in  Figure  A-2.  Since  €°  and  K  are  not  functions  of  z  but  are  mid- 

xx 

die  surface  values,  they  can  be  removed  from  under  the  summation  sign. 
Then  this  equation  can  be  written  as 


M 

x 

where 


[*]{•:}■ 


k=i 


ij /k  (zk 


N  {*„} 

2  , 

-  and 
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Figure  A-l.  Moment  on  Laminate  Beam 


Figure  A-2.  Geometry  of  an  N-layered  Laminate 


“y-Tg'Vk  ul  '  Vi> 


are  the  bending  stiffnesses  and  are  the  coupling  stiffnesses 
and  imply  coupling  between  laminate  bending  and  extension.  Therefore 
since  the  beams  are  assumed  to  undergo  only  bending,  all  =0  . 

Then  the  moment  equation  becomes 


M 


MM 


where  Du  ~  3  S  ^l^k  (zk  -  *kr-l) 


The  engineering  constants  for  the  beams  are 


llx  »  2.3  x  10  psi 


*>  1.6  x  10  psi 


G12  =  (.7-. 8)  x  10  psi 


12  *=  .3 


21  *  .01 


Then  the  reduced  stiffnesses  become 


Qll“ 


1-  V  V 
12  21 


=  2.3069  x  10  psi 


12 


*12  E2 


=  4.8144  x  10  psi 
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<22 


V  V 
L  12  21 


<>66  “  G12 


1.6048  x  10  psi 


.75  x  10  psi 


Only  one  transformed  reduced  stiffness,  Q^>  is  needed  to  calcu¬ 
late  the  stiffness.  It  is  related  to  the  reduced  stiffness  in  the 
following  manner: 


<n 


cos  6 


+  2(Q.„  +  2Q, 


,)sin^@  cos2#  + 


where  0  =  angle  from  the  x-axis  to  the  axis  the  fibers  are  directed 
along.  • 

For  the  0  degree  plies,  =  2.3069  x  102  psi  . 

For  the  —  45  degree  plies,  =  7.1592  x  10^  psi  . 

For  the  90  degree  plies,  =  1.6048  x  10^  psi  . 


Then  Dn  =  5.8233  x  102  lb-in2  . 

*2 

Since  M  -  El  =  El  K 

x  ^2 

ox 

El  K  *  D ..  K  or  El  =  D., 

X  11  X  11 

and  El  =  5.8233  x  102  lb-in2  . 


x 


then 


For  a  beam,  the  moment  of  inertia  is  I 


Using  the  typical  dimensions  b  =  1  inch  and  h  =  .085  inch  for 
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/ 


these  specimens,  I  =  5.1177  x  10-^  in^  . 

The  effective  Young's  Modulus  becomes  E  =  1.138  x  10^  psi 
for  this  composite. 


52 


Appendix  ji 
Material  Stiffness 

In  order  to  verify  the  calculation  of  stiffness  for  the  composite 
specimens,  an  attempt  was  made  to  measure  the  stiffness  experimentally. 

Weights  were  suspended  from  the  free  end  of  a  cantilevered 
specimen  and  the  tip  deflection  was  recorded.  The  stiffness  was  then 
calculated  from 


where  P  =  force  applied  at  distance,  1,  from  the  root.  Results  are 
presented  in  Table  XV. 


Table  XV 

Static  Stiffness 


P,  lb 

<5t»  in 

Measured 

El,  lb-in 

.5 

.06 

378.83 

1.0 

.12 

378.83 

1.5 

.18 

378.83 

2.0 

.24 

378.83 

The  measured  values  are  in  poor  agreement  with  the  computed  value  of 
El  ■  582.33  lb-in.  In  view  of  the  fact  that  the  computed  stiffness  led 
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to  a  predicted  resonant  frequency  of  1/  =  27.85  cps  which  was  in  much 

better  agreement  with  measured  values  (26.1  -  28.6  cps)  than  the 

value  of  V  =  20.42  cps  calculated  using  Eq  (9a)  and  the  measured 
K 

stiffness,  the  computed  El  was  accepted  as  being  correct. 
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Appendix  £ 


In  order  to  determine  how  much  of  the  total  damping  was  air 
damping,  an  aluminum  paddle,  three  inches  square,  was  attached  to 
a  specimen,  centered  and  flush  with  its  tip.  This  increased  the 
surface  area  by  a  factor  of  two.  Free  vibration  measurements  were 
made  of  the  beam  with  the  paddle  and  of  the  beam  alone.  The 
resonant  frequency  was  approximately  the  same  for  both  configura¬ 
tions. 

Air  damping  is  proportional  to  the  area  of  the  vibrating  object 
and  the  square  of  the  maximum  amplitude.  Linearizing  the  mode  shapes 
of  both  configurations  led  to  approximating  the  air  damping  in  terms 
of  an  average  amplitude.  For  the  beam  alone,  expressed  in  terms  of 
the  tip  amplitude,  this  was  0.55^.  For  the  beam  and  paddle,  the 
average  amplitude  was  0.8^.  This  yielded  a  damping  value  for  the 
beam  and  paddle  system  six  times  that  of  the  beam  alone. 

The  measured  damping  was  greater  for  the  beam  and  paddle  than  for 
the  beam.  Assuming  material  damping  was  the  same  for  both  cases,  it 
was  eliminated  from  the  difference  of  the  two  measured  values  and  the 
air  damping  was  solved  for.  The  air  damping  was  on  the  order  of  five 
percent  of  the  total  damping  of  the  beam.  Therefore,  the  damping 
results  in  this  study  are  indicative  of  material  damping  behavior. 
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